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Abstract 

A guide to the composition of the enigmatic /o(980) and ao(980) states is their for- 
mation in (^-radiative decays. Precision data are becoming available from the KLOE 
experiment at the DA<3?NE machine at Frascati, as well as results from SND and 
CMD-2 at VEPP-2M at Novosibirsk. We show how the coupling of the jo (980) to 
this channel can be extracted from these, independently of the background provided 
by a production. To do this we use the fact that the behaviour of both the /o(980) 
and a cannot be determined by these data alone, but is strongly constrained by ex- 
perimental results from other hadronic processes as required by unitarity. We find 
that the resulting coupling for the (ft - > 7/0 (980) is ~ 10 -4 GeV with a background 
that is quite unlike that assumed if unitarity is neglected. This provides an object 
lesson in how unitarity teaches us to add resonances. Not surprisingly the result is 
crucially dependent on the pole position of the /o(980), for which there are still size- 
able uncertainties. At present this leads to an uncertainty in the (ft — > /07 branching 
ratio which can only be fixed by further precision data on the /o(980). 
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1 Introduction 



Much has been written about the possible nature of the /o(980) and ao(980) states. Are 
they conventional qq scalars, or KK molecules, or 4 quark states? Are these two states with 
almost degenerate masses closely related, or is the degeneracy just an accident? 

A clue to the composition of these states is provided by the way they appear in 0-radiative 
decays. The provides us with a clean ss system, which picks out the components of each of 
these states that couples to strangeness, even though their masses are below KK threshold. 
Consequently, the /o(980), in particular, appears as a peak in the im decay distribution, as 
seen in if) — > 0(7T7t) or D s — > tt^titi), while in channels in which the / (980) is produced from 
predominantly non-strange quarks, as in ttti — > tttt scattering or in central dipion production 
with pion or proton beams, it appears as a dip or a shoulder. 

The absolute values of the branching ratio for <f> —>■ fo/aoj, as well as the ratio of these 
ratios, has been advertised in the past as a guide to distinguish nn from ss from KK and 
qqqq systems, as summarised in Table 1 PjJ |2] . 



Composition 


BR(0 -> 7 /o(980)) 


qqqq 


0(1(T 4 ) 


ss 


0(1(T 5 ) 


KK 


< 0(1(T 5 ) 



Table 1: Predictions for the absolute rate for — > 7/0 (980) depending on the 

composition of the / (980) P El- 
According to this, the first results from the two VEPP-2M experiments at Novosibirsk PHUE] 
indicate a four quark composition. However, it is now known that the strong suppression of 
the — > /q7 branching ratio in the molecular models indicated in Table 1 is wholly because 
of over-simplified modelling of the decay process in this pointed out by Oiler j^j, so 

any conclusion about a four quark composition is far from certain. 

The r]7r-channel by which the ao is identified is clean with very little background, either 
background from competing 77 final states, but also from other hadronic states — none are 
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nearby. In contrast the Jo (980) — > 7r°7r° decay gives a final state that has an experimental 
background from p°7r°, where the p° — ► 7r°7, as well ctS ct SI gnificant contribution from 
the production of a — > 7r°7r°. These backgrounds mean that the extraction of the — > 
fo (980)7 signal depends on our ability to separate the different components. The pn signal 
is characterised by a different angular dependence than that for /o7- The former is flat in 
cos8, while the latter has a (1 + cos 2 #) distribution. The KLOE data [7] reveal that the pir 
component is essentially negligible and we will take that as read. To separate the different 
scalar tttt components, in particular the a from the /q(980), models of these two states are 
used. Since the 0-radiative decay data alone cannot determine the masses and widths of 
these resonances, their key parameters are taken to be fixed by other experiments. 

For instance, the KLOE collaboration have made a determination of the BR((j) — > /07) 
from their own data and found this to be 4.5 • 10~ 4 . So why is the reanalysis presented 
here necessary? Firstly, the KLOE group adjust the / (980) mass and width to give the 
best fit to their data. This results in a mass of 969 MeV with a width of 150 — 250 MeV. 
This width is, as we shall recall shortly, far larger than other determinations of the /o(980) 
width, which are typically 50 MeV (see for example Ref. jH]). Consequently, their fo 
peak sweeps in a larger range of tttt production and this inevitably leads to a very large 
value for the BR(<p — > foj). Clearly a more conventional fo width would mean that the 
significant production of two pions with mass of 800 — 900 MeV must come from some other 
source than the /o(980). Having adjusted the parameters to maximize the agreement with 
their data, the KLOE group then add a o contribution, where the parameters for this low 
mass 7T7T enhancement is taken from the maximum likelihood fit by Fermilab-E791 to their 
D — > 3tt data. The E791 analysis is cast into doubt by the implied phase of their S"-wave 
tttt production amplitudes not being consistent with the phase determined in the tttt elastic 
scattering. Moreover, KLOE just add this contribution to that of their broad /o(980) with 
no regard for the quantum mechanical bound that probabilities must not exceed 100%, i.e. 
they violate unitarity. The KLOE data alone cannot determine the parameters of both the a 
and of the /o(980). Consequently, we embark on a new determination of BR(<fi — > /07) which 
automatically embodies coupled channel unitarity and is in accord with meson scattering 
information on txtx and KK production. We than use the radiative decay datasets to fix 
the relevant couplings. 

If the resonant states that contributed to the decay — > '-/tttt were narrow and well-separated 
then the branching ratio to any of these could be established by simply fitting data with 
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appropriate Breit-Wigner forms. However, the scalar / (980) overlaps with the notoriously 
broad /o(400 — 1200) (or a), and couples strongly to the KK channel that opens within its 
width. Moreover, since below 1 GeV the irir channel is effectively the only open hadronic 
channel, unitarity is a particularly powerful constraint. This requires all processes in which 
a dipion system is produced to have closely related final state interactions. A relationship 
embodied in Watson's famous theorem. Consequently, the S— wave nn system produced in 
0— radiative decay is already exposed in other hadronic processes, particularly high energy 
reactions where the di-pions are either peripherally or centrally produced. Only if the states 
in this system were narrow and well-separated would fitting be a matter of taking Breit- 
Wigners with couplings to be determined as the KLOE group do. 

The coupling of a resonance is only determined independently of any other resonance with 
which it may overlap, or independently of any non-resonant background, by continuing the 
resulting <p — > jim amplitude into the complex s plane to the fo pole on the appropriate 
unphysical sheet. The residue of the pole gives the coupling and this is the only model- 
independent parameter that can be determined. In the case of a narrow isolated resonance 
this coupling given by the pole residue is directly related to a branching ratio. Where 
states are not narrow, or overlap with each other or with strongly coupled thresholds, the 
branching ratio is only related to this coupling in a model-dependent way, unless one has 
data of infinite precision. Consequently the — > 77T7t branching ratios we quote are (like 
those of all authors) subject to model dependence of the fo line shape. 

2 Unitarity and underlying hadronic amplitudes 

Unitarity has to be respected. The natural way to ensure this is to represent the purely 
hadronic scattering amplitudes with definite angular momentum by the T-matrix and express 
this in terms of a real X-matrix. Then with p the diagonal phase space matrix, we have 

T — K (I — i pK)- 1 . (1) 

The amplitude F describing the production of the same final state with the same quantum 
numbers is then related to the T-matrix, provided the "production" process has no other 
strongly interacting final states, or as in an isobar picture one assumes only isobar interac- 
tions with negligible 3 or more body interactions. The relation between F and T is then 
implemented using either the P-vector or Q- vector (here called equivalently the a— coupling 
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vector), so that 

F — P {I — i pKy 1 = aT . (2) 
There are constraints on P and a, which have to be satisfied. Firstly, 

• P and a are real vectors, 

• If the .fT-matrix elements have poles, then these must also appear in the P-vector. 

• While the latter constraint is automatically built into the coupling vector representa- 
tion, care must instead be taken of the zeros of the T-matrix elements. Thus a has 
poles to remove any real zeros of the T-matrix, or its determinant. 

The P— vector and a— vector formulations are, of course, equivalent. They each embody 
the universality that demands that poles of the S'-matrix transmit to all processes with the 
same quantum numbers in exactly the same position. That this is a fundamental S'-matrix 
principle has been questioned in a series of papers by the Ishidas Pj, as being incompatible 
with quark dynamics. However the principle of universality is a consequence of causality and 
the conservation of probability for hadronic reactions totally independently of the underlying 
dynamics of the supposed constituents of hadrons. This universality is not an optional 
constraint or a matter of debate. It is a rigorous consequence of fundamental principles 
that define what is meant by the hadron spectrum. This universality of strong interactions 
applies to ^-radiative decay, since the photon-hadron interaction is electromagnetic, and the 
</> and 7r are presumed to have negligible strong interactions. 

Though the P— vector and a— vector formulations are equivalent and translatable one-to- 
one, one form may be easier to implement in practice than the other. We will use the 
coupling vector formulation. In the energy range accessed in <p radiative decay, the only 
hadronic intermediate states that can possibly enter into ttti production are tttt and KK. 
Consequently we introduce two functions a.\ and a^, which represent the coupling of the 
to 77T7T and ^KK, respectively. Provided these functions are real for tttt energies above 
threshold, then 2-channel unitarity is satisfied by construction. 

The amplitude, F, for S— wave dipion production in <fi radiative decay is thus related to the 
two hadronic scattering amplitudes, Tn for Tin — > Tin, and T12, for Tin — > KK. '1' labels 
the 7T7T channel and '2' the KK, and time reversal invariance implies that T12 = T21. Each 
amplitude and coupling function depends on s the square of the mass of the dipion system. 
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Then unitarity requires: 



ati(s) T u (s) + a 2 {s) T 12 (s) 



(3) 



The hadronic elements typically have real Adler zeros at s — s^, and the two channel 
determinant T n T 22 — T± 2 may have a real zero at s = St- These zeros do not in general 
transmit to the decay or production amplitude, being process specific. The Adler zeros of 
the hadronic amplitudes are, in the analyses we use, taken to be at the same position so 
•Sn = si 2 = s 22 = so. Moreover, the decay amplitude for (ft —> j(inr) is expected to have 
its own (process-dependent) Adler zero at s = sa- Consequently, the coupling functions are 
parametrised to incorporate this zero, while eliminating those in the underlying T-matrix 
elements and their two channel determinant, in the following straightforward way: 



where (3 is a dimensionful constant, tt = Tu(st)/Ti 2 (st) and the Pi(s) are expected to be 
represented by low order polynomials — all real. 

For the hadronic T-matrix we will use a i^-matrix analysis of meson production by Morgan 
and Pennington which updates the AMP analysis of Ref. ^U] with the constraints on near 
KK threshold production of Ref. fIJ. This is referred to as the Re VAMP analysis. The 
hadronic amplitudes embody the information we have on the structure of the /o(600) (or a) 
and the / (980). Even though this is an old analysis, the features of the embedded / (980) 
(pole position, mass and width) are in excellent agreement with the most recent data on D s 
decays into pions. The interpretation of these in terms of possible quark structure is then 
to be elucidated by the coupling of the (ft — > 7/0- 

Later we will compare the fits obtained using this set of underlying amplitudes with those 
found by using a second input, from a much more complete and much more recent analysis 
of a wide range of di- meson production data by Anisovich and Sarantsev [12]. This we 
reference as the AS analysis. The two sets of hadronic amplitudes, Re VAMP and AS have 
quite different parameters for the /o(980), see Table 2, and this difference leads to distinct 
4> — > 7/0 couplings, as we will discuss in Section 5. 
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ReVAMP 


AS 


Pole (MeV) 


£ fl = (989 - i- 22) 


E R = (1024 - i ■ 43) 


g* (MeV) 


163 


328 


9k (MeV) 


173 


398 


r w (MeV) 


25 


98 


(MeV) 


3 


36 



Table 2: Pole position, titt, KK couplings and widths relative to the / (980) 
embedded in the two different schemes ReVAMP fTI] and AS Note that in 
the AS underlying amplitudes the /o(980) is a much broader resonance than in 
the ReVAMP amplitudes. 

It is important to appreciate that hadronic data may well be fitted equally well by sev- 
eral iT-matrix parametrisations. The resulting T-matrix elements will be identical within 
experimental uncertainties and should provide the same poles of the T-matrix. However, 
depending on how the underlying i^-matrix is described in terms of poles alone, or poles 
plus some background, the poles of the .ff -matrix may be quite different in number and 
in quite different positions. This is of no physical consequence. However, it clearly does 
have consequences for those that identify the poles of the .ff -matrix with underlying bare 
states ^3]. R is important to realise that this is a modelling, that may or may not accord 
with QCD. Having a range of fT-matrix parametrisations of the same data recognises that 
we cannot yet calculate in detail the strong physics aspects of QCD. Only if one desires to 
imbue the i^-matrix elements themselves with significance is it necessary to ensure that the 
i^-matrix elements have the correct left hand cut analyticity, whilst fitting along the right 
hand cut. If the elements are just a convenient way of parametrising data along the unitarity 
cut, then approximating distant left hand cut effects by poles and a polynomial background 
is sufficiently general. 

Though the hadronic amplitudes ReVAMP and AS are the results of fits to similar datasets 
on 7T7T and KK production data, they have quite different parameters for the /o(980) as 
indicated by the pole positions, see Table 2. (The precise definition of the quantities g n , T n 
and gx^K will be explained later.) To understand the differences in these two versions of 
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Figure 1: The tttt underlying amplitudes as obtained from the ReVAMP parametrisation (solid 
line) and from the AS parametrisation (dashed line) are compared to two sets of experimental 
data. In the upper picture the modulus of the isospin zero amplitude T\\ is plotted on top of 
the results of the Ochs- Wagner analysis of the CERN-Munich data ^3]. In the lower picture the 
modulus squared of the tttt S-wave (which includes the contribution from the 1 = 2 channel |17j ) 
is plotted on top of the BNL-E852 data of Ref. On both plots the low energy data (circles) 

correspond to the BNL-E865 measurement of Kz± decays |16j . Differences between the data and 
the tttt amplitudes is attributed to non-pion exchange contributions to di-pion production, which 
are markedly different in the two cases. 
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the /o(980) we show the underlying tttt amplitude T\\ as computed in the Re VAMP and in 
the AS scheme, compared with the CERN-Munich [Hj, BNL-E852 15J and BNL-E856 PU 
data. By comparing the solid and dashed lines it is immediately apparent that the / (980) 
embedded in the AS amplitudes, which is responsible for the dip in the / = J = — > titt 
cross-section, is a much broader resonance than in the Re VAMP amplitudes. Moreover, as 
we can see in Table 2, in the AS scheme the /o(980) corresponds to a pole situated above 
KK threshold and has a mass larger than the mass of the </>, RcEr > m<p, as opposed to 
the Re VAMP 's /o- These differences will critically influence the calculation of the <p ~^ ifo 
coupling, obtained by continuing the amplitude F(s) into the complex energy plane on the 
second sheet to the / (980) pole, as will be evident in Section 5. 



3 The Fit 



Based on Eqs. (0111), our simultaneous fits to the experimental data are performed using the 
expression 

which relates the experimental decay rate to the amplitude F(s). Here y/s = E = M T7T is 
the invariant mass of the di-pion system; p(s) is the appropriate three-body phase space for 
the decay <fi — > 7(7T7t) 

7T 2 



p( S ) = ^3 ( m J- S )V S - 4m -' ( 6 ) 

and F(s) is the complex amplitude, which is related by unitarity to the hadronic amplitudes 
Tii(s) and T 12 (s) through Eq. (3). The data, of course, know about gauge invariance of the 
photon field and so the fits require the coupling functions ctj(s) to have a zero at s = m^. 
However, the fitting is aided by building in this consequence of gauge invariance from the 
start and so we introduce a factor of (1 — s/rrv^), where here is the value of the cm. 
energy to which the initial e + e~ beams are tuned (see Ref. ^S] for details). Consequently 
we construct the decay amplitude of Eq. (jSJ) using 

o.w = ^Uw + r^H 1 -^) - 

(s-so) (s-s T ) J V mil 



o» W = ^Um-P*) , (7) 

where (3 together with the coefficients of the polynomials p\{s) and P2{s) are the parameters of 
the fit. The AS hadronic amplitudes have no zero in their two channel T-matrix determinant. 
Consequently, in this case j3 = 0, or equivalently St — > oo. 
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Figure 2: Simultaneous fits to KLOE [2] and SND data as obtained by using 
the Re VAMP [HUE] set of underlying amplitudes, with constant Pi(s). 

The position of the Adler zero required by chiral dynamics is expected to be sa = 0{ m V)- 
However, our fits prove insensitive to its exact position, so we simply fix sa = 0. 

The polynomials Pi(s) are in general to be as little structured as possible consistent with the 
fitting of the data. Using the Re VAMP set of underlying hadronic amplitudes, a satisfactory 
fit is obtained by simply using Pi(s) of zeroth order in s, i.e. constant 

Pi(*)=ai, (8) 
p 2 {s) =a 2 - 

The line shape obtained from this fit, labeled by ReVAMP-0 (where indicates the order 
of the polynomial in the a's) is shown in Fig. |21 together with the data points from KLOE 
[7] and from a recent reanalysis of SND jS]. All the fits are performed by integrating the 
thoretical expression over each experimental mass bin. This is important where the ampli- 
tude and phase space are varying rapidly within a given experimental bin. The values of x 2 
and free parameters determined by this fit are shown in the first column of Table E3 
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ReVAMP-0 


ReVAMP-II 


AS-II 


X 2 /d.o.f. - KLOE 


1.60 


0.82 


1.28 


X 2 /d.o.f. - SND 


0.52 


0.62 


0.59 


ai • 10 2 


0.27 


0.51 


0.22 


61 • 10 2 




-1.19 


-3.08 


ci • 10 2 




0.83 


4.89 


a 2 • 10 2 


-0.07 


-3.44 


-10.09 


6 2 • 10 2 




11.94 


28.43 


c 2 • 10 2 




-10.45 


-14.24 


/5 • 10 2 


0.21 


0.13 





Table 3: Values of the free parameters as determined by our fits to KLOE + SND 
data corresponding to the Re VAMP and AS parametrisation for the underlying 
amplitudes Tn and T± 2 , and to a polynomials of either zeroth or second order. 



In contrast, when underlying amplitudes from the AS set are used, second order Pi(s) poly- 
nomials are needed to reach a quality comparable to that of the ReVAMP-0 fit: 

pi(s) = ai + 61 s + ci s 2 , , . 

p 2 (s) = a 2 + b 2 s + c 2 s 2 , 

as % 2 /d.o.f. is still as big as 4 for first order pi's. 

The Re VAMP amplitudes provide an excellent fit to the radiative decay data, far better 
than the comparable fits using the latest AS amplitudes. On chi-squared grounds alone 
the latter would be ruled out as quite improbable. However, we are aware that the AS 
amplitudes are the result of studying a far more comprehensive set of hadronic data than 
the older analysis included in ReVAMP and it is for this reason and that alone that we 
include the results using the AS amplitudes. Since the AS and ReVAMP differ significantly 
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Figure 3: Simultaneous fits to KLOE [Zj and SND data as obtained by 
using the ReVAMP [TT] (solid line) and the AS [12] (dashed line) underlying 
amplitudes, with Pi(s) polynomials of second order in s. 



for complex values of energy away from the real axis where our fits are performed, see Fig. 
these differences have an appreciable effect on the coupling g$ of <f> — > 7/0, as we will discuss 
in more detail in the Section 5. 

In the second and third rows of Table El we compare the values of x 2 /d.o.f. and of the 
free parameters as determined by our second order fits relative to each choice of underlying 
amplitudes, indicated as ReVAMP-II and AS-II, and in Fig. |3]the line shapes obtained from 
these fits are compared with the experimental data from KLOE [7j and SND 

Notice that global fits to all four sets of available data on <fi radiative decay, KLOE [7j, 
SND jH |3] and CMD-2 [3] are dominated by the smaller error bars from KLOE and the 
reanalysed SND , so that these fits are indistinguishable from those of Fig. |21 and El 
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4 Determination of the Couplings 



The coupling which governs the decay of 4> into /o(980), is calculated by continuing the 
amplitude F(s) into the complex s plane to the position of the fo pole (which is determined 
by the underlying amplitudes' denominator). This procedure ensures the absence of any 
background contamination. The pole residue is then evaluated from the strong interaction 
amplitudes T(s), and the radiative decay amplitude F(s), which are both pole dominated 
in the neighbourhood of the resonance pole itself, i.e. for s ~ s R 

Tii 00 = (10) 
sr — s 

F( 8 ) = (11) 
sr - s 

Thanks to the parametrisation of the T's in terms of the .ff -matrix, we know their numerators 
and denominators 

D{s) D{s) 



and consequently we have 



_ ai(s)N n (s) + a 2 (s)N K (s) 
[) ~ D(s) ' 1 ' 



In the region nearby s = sr we can write a Taylor expansion of the function D(s) truncated 
at the first order 

D{ S ~ Sr) ~ D(S R ) + D'(sr) (S - Sr) = D'{ S r) (s - Sr) , (14) 

since D(sr) = at the resonance pole. By substituting this into Eq. ()12j) and comparing 
with Eqs. JUJ and ((TTJ), we find 



9-k 



\ D'{s R ) ' (15) 

1 a 1 (s R )N 7r (s R ) + a 2 (s R )N K (s R ) 
94> = F^T — \ • 16 

9n D'(Sr) 

Knowing the — > /o7 coupling is the key result of our analysis shown in Table 0] It 
is this that is model-independent and so can be compared with predictions from different 
modellings of the make-up of the /o(980). 
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ReVAMP-0 


ReVAMP-II 


AS-II 


g<j> ■ 10 4 (GeV) 


6.21 


6.50 


18.66 


BR{(j> - / o7 ) • 10 4 


0.31 


0.34 


1.92 



Table 4: Values of the branching ratio BR(<p — ► / (980)7) an d °f the couplings 
corresponding to different fits to KLOE and SND data using the ReVAMP 
and AS ^2J underlying hadronic amplitudes. Note that BR(<f) — > /o7 ~~ ► 7r°7r°7) 
has been multiplied by a factor three to obtain BR((p — > / (980)7) to take into 
account the contribution of the charged pions decay mode of the / (980). 

4.1 Suppression effect of the photon momentum 

There is one delicate but crucial point we want to underline. What we are studying here is 
the radiative decay of the <fi. This means that the corresponding amplitude F(s) includes 
a factor given by the momentum of the radiated photon k oc (1 — s/mV), which we have 
included in the a coefficients, see Eq. (JJJ). Now, it is easy to see how this factor produces 
a very strong suppression in F(s) in the energy region around 1 GeV, exactly where one 
calculates the value of the — > / 7 coupling g^. To make it clear we construct the amplitude 
F(s) which one would obtain if the photon momentum factor was divided out of the a's. 
The plots in Fig. 0] illustrate this suppression. Let's see the details. F(s) is constructed 
from Tii(7T7r — > 7T7r) and Ti\{KK — > 7r7r) with weights given by the coupling functions 
ol\ and a 2 respectively, according to Eq. (jSJ). In the region around 1 GeV, the amplitude 
|Tn|, Fig. EJa), has a sharp dip, whereas |Ti 2 |, Fig. Efb), shows a pronounced peak (both 
structures signal the presence of the /o(980)). Because of the peak in the 1 GeV energy 
region of the (p — > /o7 spectrum, one expects the couplings a±,a2 make the contribution of 
\T 2 i(KK — > 7T7r)| dominate over |Tn(7r7r — > 7r7r)| in that region. Though F(s) seen in Fig. 4c 
shows the dominance of a^T^i, QED gauge invariance strongly suppresses the amplitude for 
small photon energy. This leads to a relative enhancement of the aiTu contribution away 
from the narrow /o(980) peak that is essential for obtaining a good fit, as seen by comparing 
Figs. 4c and 4d. 

This same photon momentum effect impacts on the determination of the value of the coupling 
and of BR(<p —>■ /07), making them both extremely sensitive to the exact position of the 
/o(980) pole, as we emphasise in the next section. 
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Figure 4: These diagrams illustrate the suppression of the <p 7 71 " 71 " decay 
amplitude, F(s), as a function of dipion mass y/s due to photon emission. Plots 
(a) and (b) show the ReVAMP hadronic amplitudes |Tn| and |T 12 | from which 
F(s) is constructed, according to Eq. (jHJ). Plot (c) shows the amplitude |.F(s)| 
one would obtain if the photon momentum factor k oc (1 — s/rrA) was divided 
out from the a\^- Plot (d) shows the amplitude |-F(s)| itself, as determined 
by our fit using the parametrisation of Eq. (J7J) for the ax,2 coefficients. It is 
immediately clear that the photon momentum dramatically changes the dipion 
mass distribution. Its effect is particularly striking in the region around 1 GeV, 
where the peak in F(s) is very narrow and spiky. Moreover, while slopes 
down gently beyond 1 GeV, plummets to zero as k — > 0. This illustrates 

why the exact position of the /o(980) pole has such a strong effect on the coupling 
and branching ratio for <fi — > 7/0 (980). 
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5 Widths and Branching Ratios 



Most authors, instead of specifying the (p — > /o7 coupling, quote the partial widths and finally 
the branching ratio BR(<p — > fol)- Then to calculate widths and branching ratios one has 
to rely on model dependent assumptions about the resonance line-shape. For instance, one 
could simply assume 

T<p = P<t>(rn 2 R ) \g<f>(s R )\ 2 (17) 

using the value of as determined by the pole residue, where the pole position is parametrised 
by Sft = m\ — imnT t . Nevertheless, to take into account the features of F(s) in the neigh- 
bourhood of the fo resonance pole (in particular its finite width), we use a more refined 
formula which expresses the width T as an average of the product of the phase space 
times the square of the coupling integrated over the f resonance density 

IX Js wn \S - S R \ 2 

where s nn is the nn threshold value of s, and for s < m 2 ^ (cf. Eq. (JJJ)) 

" 'ml-s)V~s. (19) 



2ml 



This would be the obvious procedure for integrating over the finite width of a resonance, if 
the resonance was isolated. But here the tail of the fo is not given by any Breit-Wigner shape, 
since it is intrinsically linked with the a as required by unitarity for overlapping resonances. 
Consequently, data never show a pure fo line shape. So in this case the definition of widths 
and branchings involving the fo introduce a strong model dependence in the calculation. 
Once again we stress that it would be wiser to compare values of the the — > /o7 coupling 
calculated in different schemes as the residue of the resonance pole rather than relying on 
the branching ratio. 

Notice also that whether one uses the definition of decay width from Eq. ()17j) or from Eq. (|18|) 
is optional for the fits based on Re VAMP amplitudes, but the choice of Eq. (jl8)l is inevitable 
when AS underlying amplitudes are used, since there the / pole occurs at an energy larger 
than m<£, where p^(s) = 0. A suitable normalization has to be chosen for the integral in 
Eq. ()18j) : we use the T t defined from the imaginary part of the pole position for the /o(980), 
which ensures consistency in calculating all the relevant partial widths in this process: T^, 
r,r and ry . We do not renormalize our partial widths so that the sum of + Tk so defined 
equals T t from the pole position as, for example, the authors of Ref. [12] do. 
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Table |U shows the branching ratio corresponding to the pole-determined couplings, g^, of 
the two different fits. To understand why the values differ so much, even though the fits 
look very similar, we have to take a step back and understand the differences between the 
Re VAMP and the AS underlying amplitudes. 

There the pole position relative to the /o(980) resonance is identified by continuing the T 
amplitudes into the complex plane, whereas the partial widths are obtained from formulae 
analogous to Eq. (JTSJ) 



As we discussed in Section 1, the /o(980) is a much broader and heavier object in the AS 
amplitudes than in Re VAMP 's. Now, the effect of such a broader resonance in the AS 
scheme is not that evident on the real axis, but becomes very relevant when moving into the 
complex s plane past the pole. This effect is transmitted to F(s) = a{Tu + a 2 7i2- On the 
real axis, where the fits to experimental data are performed, the F(s) amplitudes obtained 
in the two different schemes look very similar, but the similarities fade away as soon as 
we leave the real axis to move into the complex plane, where the couplings are computed. 
This is demonstrated in Fig. El where the modulus of F(s) is plotted against the real and 
imaginary part of E = */s in a 3-dimensional plot. Slices of \F(E)\ at constant values of 
ImE show the profile of \F\ when moving further and further into the second Rieman sheet, 
as the pole corresponding to /o(980) takes shape and becomes the dominant feature in that 
energy range. By comparing the two plots one can readily see the differences between the 
two cases: for Re VAMP the pole is narrow and spiky and it occurs below the mass of the </>, 
very close to KK threshold and not very far into the complex E plane. On the contrary, the 
AS pole is a much broader structure, which occurs at energies higher than and twice as 
far from the real axis. These crucial differences explain why the couplings g$ and branching 
ratios BR((p — > mrj) (shown in Table HJ) calculated with two sets of underlying amplitudes 
differ so much. 




(20) 



(21) 



where s W7r and s K -^ denote the tttt and KK threshold values of s, and 




(22) 



(23) 
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Figure 5: Three-dimensional representation of the amplitude \F\ as a function of energy, 
E. This plot shows how the amplitude changes when continued to complex values of E past 
the /o(980) pole on the second Rieman sheet. Note that in the AS scheme (lower diagram) 
the /q(980) appears as a broad and pronounced peak, which shows up over more than one 
slice. On the same scale, the Re VAMP / (980) corresponds to a much narrower and lighter 
resonance with a much more localized effect in the profile (upper diagram). 
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We can now compare the values for the g$ coupling and the — > /o7 branching ratio we find 
from our analysis to the numbers obtained by various modellings of the <fi decay into /o7- 
First we consider the decay as if it only occurred through KK loops, i.e. <fi — > KK'-f — > /07, 
as computed in the models of Ref. [TJ EH El ©• By applying the earlier models of 
Ref. PQI2], in correspondence to the / pole positions given by Re VAMP and AS, one would 
indeed find very small values of the <ft ~ > /o7 branching ratio, of order < 0(1CT 6 ), consistent 
with Table 1. In a more recent K + K~ loop calculation, Markushin presents a study of 
<p — > 7T7T7 decay in a coupled-channel model containing the tttc, KK and qq channels [20j. 
There he finds an /o(980) dynamically generated pole (i.e. a pole which disappears in the 
large N c limit) corresponding to a molecular-like KK state. Nevertheless, the branching 
ratio he quotes is BR((f) — > 7r°7r° 7) = 3.5 ■ 10 4 , of the same order of magnitude that older 
calculation would ascribe to a four-quark state ! Similarly, in Ref. jB] Oiler presents the 
latest update of a calculation of the 4> — > / 7 branching ratio [2T] based on a KK loop 
model for <ft decays into two pseudoscalar mesons unitarising Chiral Perturbation Theory. 
In this paper though [0], a <fi^K K Q contact vertex is added to the KK loop diagrams and 
finite widths effects of the intermediate resonance are taken into account. This leads him 
to find a branching ratio BR((p — > / 7) even larger than Markushin's value. However, Oiler 
notes, as we have, that the cubic dependence of the decay width on photon momentum can 
alter the branching ratio by a factor two if the resonance mass is varied by just 5 MeV! 

Our analysis is more general than these specific meson loop calculations, since it does not 
restrict <fi to decay through a particular channel, but it allows for any possible decay mech- 
anism. Nevertheless, in the limit in which only the T 2 \{KK — > ixix) underlying amplitude 
is used in Eq. (J3J) for constructing F(s), we can expect our calculation to somehow reflect 
Refs. modelling. In this particular case, we find BR{<p /07) = 0.6- 1(T 4 

for Re VAMP and BR(<p -> /07) = 1.9 ■ 1(T 4 for AS, cf. Tabled Although the fits do not 
appear too bad when plotted against the data, their x 2 /d.o.f. are in the range of 4 — 5, so 
they are not included in Table 

Finally, we compare our results to those reported in the data analyses by the KLOE [7] 
collaboration at DA$NE and by the CMD-2 pj and SND (U IE] collaborations at VEPP- 
2M. These are summarised in Table El KLOE, CMD-2 and SND approach to the problem is 
based on fitting dBR((f) — > 7r°7r°7) / dM n7T data with some appropriate Breit-Wigner forms for 
the fo and the a resonances. The branching ratio is then computed as the area underneath 
the curve determined by the fit. All their evaluations give similar results because they all 
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BR{<f> - / o7 ) • 10 4 


KLOE 


(4.47 + 0.21) 


CMD-2 


(2.90 + 0.21 + 1.54) 


SND 


(3.42 + 0.30 + 0.36) 


SND - reanalysis 


(3.5±0.3«|) 


ReVAMP-0 (fit to KLOE+SND data) 


0.31 


ReVAMP-II (fit to KLOE+SND data) 


0.34 


AS-II (fit to KLOE+SND data) 


1.92 



Table 5: Values of the branching ratio BR(<p — > /o7) as calculated by the 
KLOE [7j, CMD-2 [3| and SND jUE] collaborations compared with those obtained 
in our determination using the ReVAMP JT] and AS ^2] underlying hadronic 
amplitudes. 

closely follow the same prescription, Ref. [13] : again, decays into /o7 are modelled as if 
they proceed through K + K~ loops only. The mass spectrum is represented as a sum of three 
terms, one representing the decay of <j) i n two pions through a scalar resonance, one taking 
into account the pre background and one corresponding to the interference between them. 
The mass of the f together with the q^k+k- and <7/ o7r +7r- couplings are free parameters 
to be determined by the fit, whereas the mass and width of the a are fixed. Once again 
we stress that modelling the irir S-wave spectrum, characterized by wide overlapping and 
interfering resonances, with a Breit-Wigner model is totally inconsistent with unitarity. In 
fact, the values of the q^k+k- and gf Q7T + n - couplings they find are very large, leading to a 
disproportionately wide / (980): r w ~ 150 - 250 MeV! 
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6 Conclusions 



We have shown that there is a model independent way to determine the — ► /o7 coupling g$, 
which we believe is the only rigorously correct way to treat S'-wave interactions, where the 
underlying resonances are far from Breit-Wigner like. It is based on unitarity, which gives 
strong constraints below 1 GeV where the tttt channel is effectively the only open channel 
and it ensures universality of dipion final state interactions through Watson's theorem. The 
results for g^ are given in Table 0] and for the branching ratio in Table 03 

0-radiative decay to the / (980) in principle provides a way of determining the internal 
composition of this enigmatic scalar. Models give quite different predictions for a qqqq or 
ss state, with a range of predictions for a KK molecule. Long before experimental data 
on 0-radiative decays became available, Achasov [I] stressed that they would reflect QED 
gauge invariance with the decay distribution being proportional to the cube of the photon 
momentum. He further emphasised that models must incorporate this behaviour too [TS] . 
We have seen that these momentum factors make 0-radiative decay a very difficult tool for 
unravelling the structure of the / (980). The strong suppression as Tin mass approaches 
makes the — > 7/0 coupling, g^, very strongly dependent on the /o(980) mass and width. 
At present, analyses of the underlying hadronic processes allow sizeable variation in these 
parameters (see TableEJ). This is particularly so because of our poor knowledge of KK —>■ tttt 
scattering, as emphasised long ago, for instance in [TTj . 

Nevertheless, the underlying dynamics of the decay is clear, tttt and KK are the only 
hadronic intermediate states relevant to spin zero interactions below the 0-mass. Conse- 
quently, the decay — > ^tttt can only proceed either by the coupling ai of — > 77T7T followed 
by 7T7T interactions or the coupling a 2 for — > jKK followed by the KK system producing 
a 7T7T final state. Given the fact that the is overwhelmingly an ss system, we expect the 
coupling a?2 that picks out the KK — > tttt scattering component to dominate. That it does 
is perfectly illustrated in Fig. 4. Recall Fig. 4c shows F(s), the — > 77T7T amplitude with the 
photon momentum divided out. This is seen to look very like T12, the KK — > nir amplitude, 
particularly in the 0.95 to 1 GeV region. The tttt — > tttt component, T n , is small but not 
negligible below 900 MeV. Multiplying the amplitude F(s) by the photon momentum factor 
k oc (1 — s/mV) (required by QED gauge invariance) to get the full amplitude suppresses the 
contribution of the KK — > tttt component with its spiky peak close to 1 GeV and enhances 
lower masses where mi — > Tin interactions dominate. For instance, the photon momentum 
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enhances the region at 900 MeV by nearly a factor 3 relative to the peak at 980 MeV, and so 
in the decay distribution, where the photon momentum appears cubed, by a relative factor 
of 24. Looking at the underlying hadronic amplitudes displayed in Figs. 4a,b shows why the 
peak in the 0-radiative decay distribution in the irn channel is so much wider than the 50 
MeV of the /o(980). Thus models that neglect the (ft — » ^{int — > -nit) contribution fail to 
reproduce the decay distribution accurately. 

The KLOE experimental integrated branching ratio is BR((p — > 77r°7r°) = (1.49 ± 0.07) • 
10 -4 . An / (980) with conventional parameters of mass 989 MeV and width 44 MeV, as 
incorporated in the ReVAMP amplitudes, gives BR(cf) — > 7/0 - ► 77r°7r°) = 0.11 • 10 -4 which 
is only 10% of this total distribution. The remaining 90% is from the decay through the 
broad /o(400 — 1200) (or a). This is large wholly because of the larger phase space and 
the photon momentum factors. The coupling to this intrinsically non-strange system is 
comparable to that expected from the branching ratio of <fi — > pir. In this case, the coupling 
of the (f) — > 7/0 is much smaller than predicted in models with qqqq structure for the /o(980), 
but rather in agreement with an ss structure or within the newly extended range for a 
KK molecule [Hj. In contrast, the underlying AS amplitudes, which embed a wider and 
heavier /o, give BR(<f) — > 77r°7r°) = 0.64 • 10~ 4 which is 40% of the experimental branching 
ratio and closer to that for a qqqq composition, but the \ 2 f° r these fits are significantly 
worse (Table 3). 

Present 0-decay data do favour the conventional narrow /o(980). Nevertheless, we do need 
to fix T 11 (7T7r — > 7T7r) and T 12 (7r7r — > KK) in the 1 GeV region very accurately before we can 
reach definite conclusions about the (f> — > 7/0 coupling and its consequences for the structure 
of the /o(980). Further direct information on T n and T 12 is unlikely, so to achieve the required 
precision we need either a careful and exhaustive analysis of the _D s -decay data, which are 
presently being accumulated, or yet higher statistics results on 0-radiative decay with fine 
resolution close to 1 GeV. In either case, these data can only fix the /o(980) parameters if 
they are analysed in a way consistent with all other sources of information in harmony with 
unitarity. The necessary data are at last becoming available. How to analyse these has been 
described here. The outcome should be clear within 12 months. 
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